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Abstract 

We present the phase diagram and critical properties of a coupled XY- 
Ising model on a triangular lattice using the mean-field approximation, the 
Migdal-Kadanoff scheme of renormalization group and Monte-Carlo simula- 
tions. The topology of the phase diagram is similar for the three techniques, 
with the appearance of a phase with XY order and Ising disorder. The 
results suggest a line of transitions belonging to the 2D-Ising universality 
class in contrast with previous data indicating a new universality class. This 
model is relevant to the columnar phases of discotic liquid crystals [such as 
hexa(hexylthio)triphenylene (HHTT)] in the limit of weak intercolumn cou- 
pling. 

PACS : 61.30.Cz, 64.70.Md, 64.60.Ak, 75.10.Hk 
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I. INTRODUCTION 



Columnar phases of discotic liquid crystals are becoming a field of increasing ex- 
perimental and theoretical activity in areas traditional to magnetic systems (modulated 
phases, stability, critical phenomena, etc.) Particular attention has been given to 

hexa(hexylthio)triphenylene (known as HHTT) compounds since they were found to be ex- 
cellent photoconductors [ 0. The photoinduced charge carrier mobility is four orders of 
magnitude higher than organic materials suitable for electronic devices and opens the way 
to improved conventional and novel applications. At present, our efforts are focused on the 
understanding of the positional and orientational orderings of the columnar phases. 

The HHTT compound, first synthesized in 1984 [ ^, has been thoroughly studied by X- 
ray diffraction [ l^-^ . These measurements indicate that two intermediate columnar phases 
are present between the low temperature crystalline K phase and the high temperature 
isotropic liquid phase I. Upon cooling from the / phase, the system enters the disordered 
hexagonal Dm phase at T ~ 93°C where the disk-like molecules form columns, which are 
arranged on a triangular lattice. The Dhd phase is characterized by intracolumnar posi- 
tional and orientational short-range order and long-range order in the triangular positions 
of the columns. As the temperature is reduced even further, the HHTT compound enters 
the ordered hexagonal columnar phase {H phase) at T ~ 70°C. In this phase, the molecules 
within the columns show long-range positional and orientational order (resolution limited 
[ 0) while maintaining the triangular arrangement of columns. Inside the columns, the 
molecules of D-^ internal symmetry are equally spaced and rotated by an angle a of ap- 
proximately 45° along the column axes [Fig. |l](a)]. The columns can be either left-handed 
or right-handed. As noted by Pontes [ |^, these observations leave only three degrees of 
freedom for each column: the helicity of the column, the overall orientational phase angle of 
the column and the vertical displacement of the column. The proposed [ § structure of the 
H phase of HHTT is characterized by the formation of a -\/3 x -\/3i?30° superlattice [Pig. 
|l](b)] in the helicity pattern and vertical positions of the molecules described as follows. One 
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third of the columns are displaced by half an intracolumnar intermolecular distance in the 
z direction (along the stacking axis). This is thought to be a mechanism to relieve part of 
the frustration due to steric hindrance. The displaced columns also have opposite helicity 
to the undisplaced ones. 

Other experiments (Rayleigh scattering [ on similar compounds have served to eluci- 
date the question of the stability [ |l| of the H phase. It is currently advocated that the H 
phase of HHTT is characterized by quasi long-range order within the columns [ [l^ . Orien- 
tational order in the low temperature region of the B. phase has been explored by ground 
state calculations [ 0], and reproduces experimental observations [|^. 

A model has been proposed recently to reproduce the sequence of phases of the HHTT 
compound [|12[- This simple model Hamiltonian incorporates, to a certain extent, what was 
noted by Pontes [ Q and mentionned earlier. It is a generalization of a model introduced by 



Plumer et al [ The columns are taken as rigid helices inscribed in a two-dimensional 
triangular array. The state of a column is determined by its helicity Ki (Ising-like variable) 
and a global phase angle 6i {XY-like variable). The vertical displacement of a column is 
then viewed as a change in the global phase angle Of. the displacement of a column in the z 
direction by half an intracolumnar intermolecular spacing is replaced by the addition of the 



angle a/2 to 6i [[T^- This simplifying approximation neglects fluctuations in positional and 
orientational degrees of freedom of individual molecules and mimicks the change in helicity 
of a column, which would result from a gradual unwinding of the helix, by a simple flip in 
the variable K^. This model is thought to be valid in the strong intracolumn coupling limit 



of the exact interaction energy [ |TT 



To second order in the moments of the mass density of the columns, the interaction 



energy of the system, coupling XY and Ising variables on a triangular lattice, is [ 112 



-?E(l-^^^.)cos(0i + 0,) (1) 

^ {ij) 
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with Ki = ±1 and (pi = 39i (0 < 0i < 27i). (We use primed parameters to be consistent with 



earher calculations [ |T2| and to leave unprimed parameters for the exact interaction energy 
[ |11[.) Note that the sign of G' is not relevant since changing G' —G' and 0j — > + 7r/2 
leaves (|I|) invariant. This Hamiltonian with G" = has been studied extensively in relation to 
the fully frustrated XY models and Josephson junction arrays in a transverse magnetic field 
[ p!^ , [T5| . They reported a single continuous transition from an XY and Ising ordered phase 
to a totally disordered phase with exponents u = 0.83(4) and P/i^ = 0.14(2) associated with 
the Ising order parameter [ |16| . These results suggested the possibility of a new universality 
class. We note also that a similar model with no Ising variables [Hamiltonian (|l]) with 
Ki = Kj = 0] has been studied recently [ p!7| . 

It is also of interest to summarize briefly the intrinsic properties of Hamiltonian (P in 



the low temperature portion of the phase diagram [ |T^. For G' = 0, the effective exchange 
parameter for KiKj brings true long-range ferromagnetic order in the Ising variables, while 
the XY variables should exhibit Kosterlitz-Thouless order. For J' = 0, the effective exchange 
parameter for KiKj is antiferromagnetic. On a triangular lattice, this results in an Ising 
disordered state at all temperatures [|T3]. On the other hand, the term cos(0j+0j) breaks the 
rotation invariance permitting true long-range order in the XY variables. The competition 
between the J'- and G"-terms are responsible for a rich T — J' phase diagram within mean- 
field theory [0. 

In this work, we study the T — G' phase diagram and critical behavior of Hamiltonian 
(H) for J' = 1. In Section ||, the mean-field phase diagram (following [ 0) is presented 
to give an indication of the type of order which may occur. In Section pT| , we present the 
Midgal-Kadanoff renormalization group study of Hamiltonian (|l]). Then, the phase diagram 
is calculated along with the critical exponents from preliminary Monte-Carlo simulations 
(Section |IV|). Finally, general remarks and conclusions are drawn in Section [V]. 
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II. MEAN-FIELD PHASE DIAGRAM 



The mean-field phase diagram for G' = —1 has already been published [ 12|. Here, we 
use the same technique to calculate the J' = 1 phase diagram (J' = 1 fixes the energy scale). 
The main purpose of this calculation is to gain knowledge of the order parameters describing 
each phase. To sixth order in a Landau-type expansion of the free energy for a unit cell (see 
Ref. [ [T^)) the resulting phase diagram is given Fig. 

Despite the fact that the mean-field approximation is crude, a lot of information can be 
extracted from Fig. ^. The low temperature phase (phase III) is characterized by XY and 
Ising order [{(pi) = {(pj) = n/2 and {{Ki) {K2) {K^)) = (+ + +)]• Phase II shows the existence 
of Xy order with Ising disorder [{(pi) = {(pj) ^ and {{Ki) {K2) {K^)) = (000)]. Finally, the 
high temperature phase (phase I) has both XY and Ising disorder. As expected, for G' = 0, 



we have a direct transition from phase III to phase I [ |Tj] , but the order of the transition is 
wrong (see below) within the mean-field approximation. The G'-term is responsible for the 
stabilization of phase II with XY order and Ising disorder. In fact, this phase II appears only 
for G' ^ 0.25. Note that italic letters (/, Dm, and H) refer to phases reported experimentally 
in HHTT and roman numerals (I, II and III) denotes of the phases in the present study. 

III. MIGDAL-KADANOFF RENOMALIZATION GROUP APPROACH 

In order to go beyond the Landau-type mean-field theory of Section || and include 
fluctuations of the order parameters, we apply the Migdal-Kadanoff renormalization group 
scheme to a wider class of Hamiltonians 

3 J' 

pnRG = E(i + cos(</). - (p,) 

^ {ij) 

-^Y.(l-K,K^)cos{(P, + (Pj) 



^^'T.K^K, (2) 
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of which Hamiltonian (|l]) is the L' = hmiting case. The term —^^(ij) KiKj is added 
to (HD since its form is not preserved under renormahzation (the renormahzation process 
creates KiKj couphngs). The Migdal-Kadanoff approach has been apphed to Hamihonian 
(0) with G" = by a number of authors [ |l^,|19[. Those calculations were carried out on a 
square lattice with rescaling factors 6 = 2 and 6 = 3 in order to detect an antiferromagnetic 
Ising order at large and negative L' . The addition of the G"-term on a triangular lattice, 
representing an intrinsic Ising frustration, is responsible for the apparition of a new phase 
characterized by XY order with Ising disorder, and should be given special attention. For 
these reasons, we will calculate the phase diagram of (0), for G' < J', using the Migdal- 
Kadanoff approximation implemented by Lee et al [ |l4l, directly on the triangular lattice 
with a rescaling factor 6 = 2. 

To apply the Migdal-Kadanoff transformation, bonds are moved in such a way that half 
the sites (to be integrated out) are linked to only two neighbors (in one spatial direction). 
This bond moving leads to a one- dimensional decimation to obtain the effective exchange 
parameters between remaining sites. In terms of new variables u{9) = exp[U{6) — U{0)], 
v{e) = exp[V{e) - V{0)] and z = exp[f/(0) - 1^(0) + (3L'] with U{e) = /3J'cos(^) and 
V{9) = pC cos{9), the recursion relations are 

Z^AA<j),-<f)2)+Z-^A2 



^^^^^-^^^^ = .M,(0) + z-M,(0) 

^(01 + 02) +^4(01 + 02) 



K(01 + 02)] 



Z 



^3(0) + ^4(0) 
,,2 _ zMi(0)^-M2(0) 
A3(0)+A4(0) 



(3) 



where 



A^i9) = 


I 2^"(^-^ 


i)u\9) 




A^iO) = 


Jo ^"^^-^ 


b)v\e) 




Asie) = 


Jo 2^""^^-' 


i)v\e) 




A4{9) = 


lo 2^""^^'' 
6 


b)u\e). 


(4) 



The primed variables in are the renormalized exchange parameters (after decimation) 
and should not be confused with J', G' and L' which are the coupling constants of the 
original Hamiltonian (^. In deriving (|^), the only assumption on the potentials U{6) and 
V{9), during the renormalization process, is their periodicities. 

The phase diagram obtained by numerically iterating (^ is presented in Fig. ^ for J' = 1 
and L' = 0. It is comprised of four phases: I- a high temperature disordered phase, II- a 
phase with XY order with Ising disorder. III- a low temperature fully ordered phase and IV- 
a phase with Ising order with XY disorder. The low temperature phase III is characterized 
by XY and Ising order [f3J' (/3J')^ f3G' {f3G'y with {(3J'y = (/3G")^ and (3L' oo]. 
The Hamiltonian with f3G' = [ |14[ a is stable limit under the renormalization group 
transformations (0), as well as the f3J' = limit. Thus, in the f3G' = case, we get 
PJ' {pj'y, [3G' — > and I3L' oo. The intermediate temperature phase II shows XY 
order with Ising disorder {(3 J' {(3J'y, (3G' {(3G'y with {(3J'y = {(3G'y, and (3L' 0). 
The PF line (Fig. ^) corresponds to an Ising-like transition within this approximation. In 
fact, putting {(3J'y = {f3G'y in (g) and (|), we find A*{0) = A*{0) and 

l^f = '-Af- (5) 

which is the Migdal-Kadanoff recursion relation for the spin-1/2 Ising model with 6 = 2. 
The CPE line is relatively difficult to locate due to the high temperature drift suffered 
by the Migdal-Kadanoff approximation. The other intermediate temperature phase IV is 
characterized by Ising order with XY disorder {(3 J' 0, (3G' — > and (3L' oo). This 
phase is not seen in the mean-field approach, and we argue that it is an artefact of the 
renormalization group scheme. In fact, by looking at Hamiltonian (g) [or with L' = 0), 
one is easily convinced that there cannot exist Ising order without XY order. The PD 
line cannot be located precisely since the Migdal-Kadanoff approach does not reproduce a 
line of true fixed points for the XY model. Nevertheless, this scheme generates a line of 
almost fixed points which gives a hint of Xy-like ordered phases. Finally, because of these 
shortcomings of the Migdal-Kadanoff approximation, the exact topology of the bifurcation 
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point P (insets A and B of Fig. ^ cannot be resolved. 

As can be observed in Fig. every phase transition is second order. The technique used 
to arrive at this conclusion is a generalization of the procedure known to reproduce the first 



order transition of the g-state Potts model for g > 4 [ Our specific implementation is 
due to Lee et al [ . The idea is to have vacancies to prevent overestimation of the local 
order in the system. 

The lattice-gas version of @) reads 

BJ' 

PHlg = -V E(l + KiKj)Utj cos(0i - <Pj) 
— ^(1 - KiKj)titj cos(0j + (pj) 

^ (ij) 

-f E [L'K^K.Ut, + K'Ut,) + AJ2U, (6) 

{ij) « 

where the vacancy variable is = 0, 1 and A is the fugacity. Applying the same renormal- 
ization scheme to (H), with w = exp[i^' + V^(0) + U{0)] and y = exp(A), we get the following 
recursion relations 

2 + - 02) + ^-'^2(01 - 02)] w^y~^/' 



[U'ih - 02)] 
K(</>l + 02)] 



2 + [zMi(O) + ^-^2(0)] w;4^-2/3 
2 + [^3(01 + 02) + ^4(01 + 02)] w'y-^/' 



\w 



2+[A3(0)+A4(0)]w7V'/' 

2 wY^^ {2 + [^Mi(O) + 2-^2(0)] {2 + [^3(0) + ^4(0)] w^^/''} 



2 ^ 2+[z'AM + z-'A2mw'y-'l^ 
^""^ 2 + [A3(0) + A4(0)]wV'/=^ 



with 



B,= r'^^v\4>). (8) 
JO 27r 

The Atj-term has been distributed equally between bonds at site i before moving them. This 
is necessary to preserve the site density at each step of the renormalization process. The line 
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of discontinuity fixed points, signaling the first order transition, appears for w,y, z ^ oo [ 
p^ , pl| . In this hmiting case, an almost fixed line appears when 

The criterion to have a discontinuity fixed point (corresponding to a first order transition) 
states that eigenvalue A, for small deviations from the fixed point along the unstable di- 
rection, must be given hj X = b'^ {d is the dimensionxnality of the lattice) [ We can 
see, from @, that this is certainly not the case. Therefore, within the Midgal-Kadanoff 
approximation and using the procedure of Nienhuis et al [ |2l| , the transitions in Fig. ^ are 
continuous. 



It is of interest compare our results with Lee et al [ |TJ] in the limit of G" = and L' = 0. 
Within the Migdal-Kadanoff approximation, both calculations give the sequence of phases 
III ^ IV ^ I with second-order transitions between them. 



IV. MONTE-CARLO SIMULATIONS 



In this section we report the results of a preliminary Monte-Carlo study of the phase 
diagram and critical behaviour of the model (|l|). We used the traditional Monte-Carlo tech- 
nique to locate the phase boundaries and to verify the overall temperature behaviour of 
the thermodynamic quantities, especially the order parameters and susceptibilities. Then, 
we used the Monte-Carlo histogram technique to locate accurately the transition tempera- 
tures and to estimate critical exponents by finite-size scaling analysis. Both Monte-Carlo 
approaches used a sequential version of the Metropolis algorithm to test separately the Ising 
and XY trial configurations. A Monte-Carlo step (MCS) is as follows. For each lattice site, 
the Metropolis test is apphed to the Ising variable {Ki) and then to the XY variable {(pi). 
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A. Monte-Carlo phase diagram 



Fig. ^ shows the resulting Monte-Carlo phase diagram. It has been obtained for a. L x L 
lattice (L = 36) with 10^ MCS for thermalization and 4 x 10^ MCS for averaging. We use 
the same convention, as in the mean-field (Section y) and renormalization group (Section 
pn] ) phase diagrams, to name the phases (roman numerals). To describe the order present 
in each phase, we define the following order parameters 

1 



K 



Mj 



L2 



(10) 

(11) 



where Sj denotes one of the three sublattices and Si = (cos(0j), sin(0j)) are pseudo-spin 
variables. The three XF-like Pmj order parameters are very sensitive to sublattice switching. 
We therefore performed the simulations with the following XY order parameter, 

Pm = max(PMi) (12) 

which has smoother temperature dependence than Pmj- In order to understand each phase 
in Fig. ^, we present the order parameters Pk and Pm as functions of temperature [Fig. 
|5|(a)]. Phase III {Pm and Pk 7^ 0) is the phase with XY and Ising order. The intermediate 
phase II keeps the XY order {Pm 7^ 0) but looses the Ising order {P^ = 0). At high 
temperatures, the isotropic I phase has XY and Ising disorder {Pm and Pk = 0). A similar 
behaviour to that of the XY order parameter Pm and its susceptibility xm [Fig- 11(b)] in 
the phase II has been reported for the six-state clock model [ ^ and is characteristic of 
systems with intrinsic disorder. 

The general features of the Monte-Carlo phase diagram can be understood easily. For 



G' = 0, there is a single transition from phase III directly to phase I [ This can be 
seen by noting that the mechanism required to have a transition in the Ising variables is 
the appearance of domain walls. For G' = 0, at such a domain wall, the XY variables are 
decoupled since (1 + KiKj) = 0. In this case, Ising disorder induces XY disorder. On the 
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other hand, for G' ^ J', there should be no Ising order in the system (Section 1). This is 
seen in Fig. ^: the PQi hne of transitions goes toward T = as G" gets larger. Ultimately, 
for G" ^ oo, there is only one transition from a phase with XY order and Ising disorder to 
the phase I. 



B. Finite-size scaling 

The rest of this Section is concerned with the finite-size scaling at the four points (M, 
N, Qi and Q2) shown in Fig. |[ The goal here is to extract only rough estimates of the 
critical exponents. We are also interested in detecting the order of the transitions and the 
possible presence of a Kosterlitz-Thouless transition. Finite-size scaling using the Monte- 



Carlo histogram method has been recently reviewed by Plumer et al [ In our case, it 



involves the specific heat (C^,), the energy cumulant {Ue), the order parameters (Pm and 



Pft-), their susceptibilities (xm and Xk)^ their logarithmic derivatives (Vm and Vk) \ ^ 
and their fourth-order cumulants {Um and Uk) [^^- We performed a simulation for each 
lattice size L = 54, 63, 72, 81, 90 and 108 at our best estimate of the transition temperature 
To from the results for the order parameters, such as in Fig. |[ Typically, 2 — 3 x 10^ MCS 
were used for thermalization and 6 — 10 x 10^ MCS were kept for averaging. Since we did 
only one simulation per lattice sizes at T^, the errors bars on the exponents are likely quite 
large (typically 10 - 15 %). 

Let us focus on the transition M, at G' = where To = 2.175 in Fig. ^, to explain 
the details of the scaling analysis used here and because these results will be compared 
with earlier published data [ |1^. The well-studied fully frustrated XY model is also closely 
connected with Hamiltonian with G" = [ p2| . 



To determine the critical temperature Tc, we used the scaling with (which assumes 
1/ ~ 1) of the positions of the extrema of the susceptibilities {xm and xk) and of the 
logarithmic derivatives {Vm and Vk)- In the special case of G' = 0, the cumulant-crossing 
method (using Uk) was also applied to get a better estimate of T^ [ U^. The results are 
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summarized in Fig. p. Both techniques give roughly the same Tc. The scahng of the 
extrema of xai and Vm gives Tc-m = 2.1631(10), while the scaling for xk and Vk gives 
Tc-K = 2.1635(7). Using the cumulant-crossing we get Tc-cc = 2.1644(10). From our data 
we can conclude that, given the errors on the critical temperatures, the transitions involving 
the XF-type and Ising-type order parameters occur at the same temperature, as expected 



22| , p!5| . The precision is quite good, especially when it is compared to the estimation of Tc 



in other systems [ . The value of Tc is required to find (3 /p because Pk ~ L'l^''' at Te. 

The other exponents can also be found by a scaling at Tc, but in our case the scaling of the 
extrema of the thermodynamical quantities proved to be more precise. 

The size dependence of the minimum of the energy cumulant (f/g) is given in Fig. |^ for 
the four transitions. The fact that t/g ~^ 2/3 when L ^ oo indicates that the four transitions 



are continuous [|30|. This is in aggreement with our renormalization group results (Section 



np and with the results of Granato et al [ |l^ in the case of G' = 0. 



Fig. I and || show the scaling that leads to the critical exponents uk and Pk/^k as- 
sociated with the Ising order parameter Pk- In the case of Pk/^^k, we used the critical 
temperature given by the scaling of the extrema of xk, Vk, Xm and Vm- The rough esti- 
mates are = 1-1(1) and Pk/^k = 0.18(1). The errors on the exponents are estimated 
only by the goodness of the fit on the plots and do not account for the (unknown) statistical 
error for each run. From the relatively large scatter of the points on those figures, it appears 
that there are large statistical fluctuations. Our estimates of the exponents are in apparent 
disagreement with those of Granato et al [ |15[. Their study involved very long Monte-Carlo 
runs but on relatively small lattice sizes {L < 40). Scaling with only our three smallest 
lattice sizes (L = 54, 63 and 72), gives uk — 0.74 and Pk/^k — 0.10, which are closer to 



their results. It is clear that the are significant finite-size effects. Fig. ^ presents the scaling 
of XKmax with L, glvlug Jk/^k = 1.55(16). 

These exponents {uk, Pk/^k and 'Jk/^k) are close to the 2D-Ising universality class 
{uk = 1, Pk/j^k = 0.125 and 'Jk/^k = 1-75). In this case, the specific heat diverges like the 
logarithm of L- Fig. ^ shows C^max against In L- The straight line fit for L > 63 is a strong 
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indication that, contrary to what was suspected earher [ |T5|, the transition at G" = is of 
the 2D-Ising universahty class. Such logarithmic divergence of Cvmax for the infinite lattice 
is also reported for the fully frustrated XY model [ . Extensive Monte-Carlo simulations 
involving large lattice sizes and better statistics are required to verify this conclusion. 

Table | lists the exponents for the four transitions considered, along with and T^. The 
missing values for the exponents at points Qi and Q2 reflects the fact that for G' = 0.4, 
the transitions involve only Ising and XY order respectively. The errors on the exponents 
are estimated only by the goodness of the fit. The exponent a = (log) means that 
Cvmax ~ In(-^v), as in the 2D Ising model. We note that 'Jm/^m at G' = 0.4 was estimated 
at Tc, because xm shows only a broad maximum across the phase II and not a sharp peak 
[see Fig. |5|(b)]. It was impossible to get an estimate of Pk/^k for transition Qi, the data 
being too erratic. 

Under the assumption that for G' = 0.0 (and G' = 0.1) there is only a single transition, 
there can be only one caracteristic length that diverges. The exponents and um are related 
to the divergence at Tc of the Ising type and XY type correlation lengths, respectively. Thus 
i>K has to be equal to um, if there is a direct transition from phase III to phase I. This seems 



to be the case given the accuracy of our results. As noted in Ref. [ 0], however, the Ising 
and XY transitions may be decoupled if z/ = 1. This question has also been addressed in 
relation to the antiferromagnetic XY model on a stacked triangular lattice [ ^ . 

Except for the values of 13x1 vk at points N and Qi, from the data given in Table | and 
Fig. 1^, we tentatively conclude that the transitions along the MPQi line are continuous 
and belong to the 2d-Ising universality class. We exclude the possibility of a transition 
of the Kosterlitz-Thouless type since at points M, N and Qi (Fig.|D, vk — — 1 and 



Cvmax is proportional to ln(L) at T^. [ |3T| , |32| . This in agreement with our renormalization 
group treatment of model (|1|) along the PF line of Fig. ^ in Section As for the other 
exponents, 'Jk/^k and 'jM/t^M seem to be stable along the MPQi line and suggest the 2D- 
Ising universality class. However, the data for Pk/^k and Pm/^m are very noisy (see Fig. 
H) and no definite conclusions can be drawn. 
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The last transition studied (Q2) is caracterized by the loss of XY order while the Ising 
variables are disordered. The specific heat scaling shows no dependence on lattice size 
indicating that a ~ 0. This rules out the possibility of a Kosterlitz-Thouless transition [ 



33| . A similar behaviour has been observed in frustrated antiferromagnetic XY models in 
a magnetic field [ In this case the exponents were reported to be nonuniversal. This 
might be the case here, but several points along the PQ2 would have to be done (with much 
better statistics) in order to follow the evolution of the critical exponents. 



V. CONCLUSION 

We have presented the phase diagram of a model for columnar phases of HHTT in the 
limit of weakly interacting columns. Its Hamiltonian couples XY and Ising variables. We 
used the mean-field approximation, the Migdal-Kadanoff scheme of renormalization group 
and Monte-Carlo simulations. These techniques give qualitatively the same topology for the 
phase diagram, showing a phase with XY order and Ising disorder. Renormalization group 
results and Monte-Carlo simulations suggest that all that transitions are continuous. From 
the renormalization group analysis, we found that the transition from phase III to phase II 
should be of the Ising type (Fig. Preliminary Monte-Carlo simulations on large lattice 
sizes confirm this result and that the transition from phase III directly to phase I is also of 
the 2D-Ising universality class. (We emphasize, however, that our statistical errors appear 
to be quite large.) This latter conclusion is in disagreement with the suggestion by Granato 
et al [ IlH] of a new universality class. 
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FIGURES 

FIG. 1. (a) Schematic side-view of a half period of HHTT in the H phase. The full period 
is composed of 8 molecules of D3 internal symmetry symbolized by triangles. The two helicities 
Ki = zizl are represented, (b) Structure of the H phase in the basal plane showing the unit cell 
(dashed line). Displaced columns are represented by open dots and undisplaced columns by full 
dots. 

FIG. 2. Mean-field phase diagram of Hamiltonian (1) for J' = 1. Different phases are noted 
by roman numerals. Phase III is caracterized by XY and Ising order, phase II by XY order with 
Ising disorder and phase I by XY and Ising disorder. Helicity and angular configurations are 
represented, and explained in the text. Solid and dashed lines represent first- and second-order 
phase transitions, respectively. The square denotes the location of critical end point. 

FIG. 3. Topological features of the phase diagram of Hamiltonian (2) within the Mig- 
dal-Kadanoff renormalization scheme for J' = 1 and L' = 0. Phases are labelled as in Fig. 2, 
with phase IV representing XY disorder with Ising order. All transitions are continuous (dashed 
lines). The structure at point P cannot be determined accurately; the two possibilities are repre- 
sented by the insets A and B. 

FIG. 4. Monte-Carlo phase diagram, with J' = 1, for a system size of 36^ sites. Different 
phases are reprensented as in Fig. 2. Full circles denote the transition temperatures for a given 
G' . All transitions are continuous. Full squares and capital letters refer to the transitions where 
finite-size scaling was performed. 

FIG. 5. (a) Temperature dependence of the Ising (Pr) and XY (Pm) order parameters for 
G' = 0.4. (b) Temperature dependence of the susceptibilities associated with Pk and Pm- 

FIG. 6. Estimation of Tc for G' = 0. The solid lines shows the scaling of the extrema of xk, 
^K, Xm and Vm vs L^^. The dashed lines represents the scaling of Tc with respect to ln~^ (L' / L) 
using the cumulant-crossing method. 
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FIG. 7. Scaling of the minimum of the energy cumulant for the four transitions considered. 
Sohd hues are guides to the eye. 

FIG. 8. Scahng of the logarithmic derivative of Pk- The solid line through the five largest 
lattice sizes gives the exponent uk = 1-1(1) • The dotted line is for the scaling using the two 
smallest lattice sizes (see text). 

FIG. 9. Scaling of the order parameter Pk at Tc = 2.1635(7). The solid line through the five 
largest lattice sizes gives the exponent Pk/i'k = 0.18(1). The dotted line is for the scaling using 
the three smallest lattice sizes (see text). 

FIG. 10. Scaling of the susceptibility associated to Pk- The solid line through the five largest 
lattice sizes gives the exponent jk/i^k = 1.55(16). 

FIG. 11. Specific heat as a function of lattice size ln(L). 
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TABLES 

TABLE I. Summary of the critical exponents for the four transitions considered. The errors are 

estimated by tlic robustness of tlic fit of the data and can be larger due to statistical fluctuations. 
The transitions Qi and Q2 involve XY and Ising orders, respectively, and only those exponents 
are given. Question marks indicate large scatter in the data. 





G' = 0.0 (M) 


G' = 0.1 (N) 


G' = 0.4 (Oi) 


G' = 0.4 (O9) 


Tn 
-*- 


2.175 


2.165 


1.745 


2.37 


To 


2.1635(7) 


2.148(2) 


1.737(3) 


2.376(5) 




1.1(1) 


1.18(16) 


1.06(12) 






0.18(1) 


0.57(1)? 






IkI^k 


1.55(16) 


1.58(16) 


1.64(16) 






0.90(8) 


1.19(13) 




1.09(7) 




0.44(5)? 


0.56(10)? 




0.58(14)? 




1.46(6) 


1.46(18) 




1.10(14) 


a 


(log) 


(log) 


(log) 
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